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SUMMARY 


The  effect  of  laminate  thickness  on  the  interlaminar  stresses  tn 
rectangular  quasi-isotropic  laminates  under  uniforir.  axial  strain  was 
studied.  Laminates  from  8-ply  to  infinitely  thick  were  analyzed.  Thick 
laminates  were  synthesized  by  stacking  (45/0/-45/90)  ply  groups,  rather 
than  grouping  like  plies.  Laminates  with  and  without  delaminations  were 

studied. 

in  laminates  without  delaminations,  the  free-edge  interlaminar  normal 
stress  distribution  in  the  outer  ply  groups  was  insensitive  to  total 
laminate  thickness.  The  interlaminar  normal  stress  distribution  for  the 
interior  ply  groups  was  nearly  the  same  as  for  an  infinitely  thick 
laminate.  In  contrast,  the  free-edge  interlaminar  shear  stress  distri¬ 
bution  was  nearly  the  same  for  inner  and  outer  ply  groups  and  was 
insensitive  to  laminate  thickness.,  In  laminates  with  delaminations,  those 
delaminations  near  the  top  and  bottom  surfaces  of  a  thick  laminate  have 
much  larger  total  strain-energy-release  rates  (G.j,)  and  mode  I-to-total 
(Gj/G.j,)  ratios  than  delaminations  deep  in  the  interior.  Therefore, 
delaminations  can  be  expected  to  grow  more  easily  near  the  surfaces  of  a 
laminate  than  in  the  interior.  This  is  consistent  with  experimental 
results  reported  in  the  literature.  Also,  near  surface  delaminations  in 
thin  laminates  tend  to  have  larger  strain-energy  release  rates  than 
corresponding  near  surface  delaminations  in  thick  laminates. 

INTRODUCTION 

interlaminar  stresses  develop  at  free  edges  and  delamination  fronts 
in  composite  laminates  loaded  in-plane  because  the  individual  anisotropic 
plies  have  different  mechanical  and  thermal  properties.  In  particular, 
differences  in  Poisson's  ratios,  shear-extension  coupling  coefficients. 


and  thermal  expansion  coefficients  cause  each  ply  to  try  to  deform  differ¬ 
ently.  Compatibility  of  ply  deformations  is  enforced  by  interlaminar 
stresses  near  free  edges  and  delamination  fronts. 

Interlaminar  stresses  have  been  studied  extensively  for  thin  lami¬ 
nates,  e.g.  4-  and  8-ply  laminates.  But  most  practical  laminates, 
especially  those  for  compression  structures,  consist  of  considerably  more 
than  eight  plies.  A  few  studies  of  thicker  laminates  (refs.  1-2)  have 
found  that  thicker  laminates  are  more  prone  to  delaminate  than  thinner 
laminates  of  the  same  stacking  sequence;  that  is,  a  s 

laminate  will  delaminate  at  a  lower  strain  level  than  a  [0/45/-45/90] ^ 
laminate.  But  like  plies  need  not  be  grouped  together  in  thick  laminates. 
If  like  plies  a^e  not  grouped  together,  are  thicker  laminates  more  or  less 
prone  to  delaminate  than  thinner  laminates?  References  3  and  4  report 
results  for  related  laminates  in  which  the  thinner  laminate  delaminates, 
but  not  the  thicker  one.  The  objective  of  this  paper  is  to  analytically 
determine  the  effect  of  laminate  thickness  on  interlaminar  stresses  in 
rectangular  quasi-isotropic  laminates  of  different  thicknesses  under 
mechanical  loading.  Laminate  thicknesses  from  8-ply  to  infinitely  thick 
were  considered.  Laminate  thickness  was  varied  by  changing  the  number  of 
plies  instead  of  changing  the  ply  thickness.  For  laminates  without 
delamination,  the  free-edge  interlaminar  stresses  were  examined  through 
the  laminate  thickness  and  across  the  laminate  width.  Quasi- three- 
dimensional  (Q3D)  finite  element  analysis  was  used  to  calculate  stresses. 
Interlaminar  stresses  near  delamination  fronts  in  delaminated  specimens 
were  characterized  in  terras  of  strain-energy-release  rates.  Strain- 
energy-release  rates  were  calculated  two  ways:  with  a  Q3D  finite  element 
analysis  and  with  a  classical  laminate  theory  (CLT)  analysis.  The  results 
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of  the  free-edge  stress  analysis  and  the  strain-energy- re lease  rate 
analyses  were  used  to  discuss  a  few  relevant  experimental  results  from  the 
literature • 

NOMENCLATURE 

a  delamination  length,  m 

b  serai-width  of  the  rectangular  laminate,  ra 

Eii  Young's  modulus  for  orthotropic  material  in  the  i-direction,  MPa 

shear  modulus  for  orthotropic  material,  MPa 
G,j,  total  strain  energy  release  rate,  J/m^ 

If  Ilf  III  strain  energy  release  rates,  respectively,  J/m^ 
h  ply  thickness,  m 

N^,Ny,Nj^y  stress  resultants,  MPa-m 

MxrMy,M^y  moment  resultants,  Nm/m 

t  total  thickness  of  the  laminate,  m 

U,V,W  displacement  functions,  m 

u,v,w  displacements  in  the  x-,  y-,  and  z-directions,  respectively,  m 

x,y,z  Cartesian  coordinates,  m 


K  ,<  ,IC 

x"  y'  xy 


uniform  axial  strain  in  the  x-direction  (e^  =  0,001) 
mid-plane  curvatures, 


''ij 

Poisson's  ratio 

M 

Cartesian  stresses. 

0  ,  a 

X 

Subscrip 

ts 

if  j 

1f2,3' 

1f2,3 

longitudinal,  transverse, 
of  a  unidirectional  ply 

I, II, III 

Mode  I,  Mode  II,  and 

Mode 

i 


others 


^ '  *  *  ■ 
i  i  ■/ 


f 

i'  V'- 


IT/, 


[(45/0/-45/90)J^ 


(45/0/-45/90/90/-45/0/45)  is  a  repaating  unit 
repeated  n  times 


[(45/0/-45/90)J^ 


[ (45/0/-45/90)/(45/0/-45/90)/...  n  timesj^ 
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DESCRIPTION  OF  INTERLAMINAR  STRESS  PROBLEMS 
The  two  basic  problems  studied  were  the  free-edge  stress  problem 
(Fig.  1(a))  and  the  edge  delamination  problem  (Fig.  1(b)).  For  both  cases 
the  laminate  is  assumed  to  be  long  and  rectangular.  Away  from  the  ends 
where  the  loads  are  applied,  the  displacements  at  any  x  «  constant  plane 
are  assumed  to  be  given  by  (refs.  5,  6) 

u(x,y,z)  =  X  +  U(yi,z) 

v(x,y,z)  =  V(y.,z)  (1) 

w(x,y,z)  =  W(y^z) 

The  €  is  the  uniform  axial  strain  and  U,  V,  and  W  are  functions 
o 

of  y  and  z  only.  Equations  (1)  describe  a  **quasi-three-dimensional" 
(Q3D)  problem.  The  modifier  "quasi"  is  used  because  there  are  displace¬ 
ments  in  three  directions,  but  the  gradients  of  U,  V,  and  W  with 
respect  to  the  x-coordinate  are  zero. 

For  symmetric  laminates  without  de laminations,  eqs.  (1)  rigorously 
define  the  displacements.  But  for  unsymmetric  laminates  or  symmetric 
laminates  with  de laminations,  there  is  usually  a  tendency  to  twist  about 
the  X  axis.  This  twist  is 

K  =  -  (2) 

xy  3x3y 
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Since  w  is  not  a  function  of  x  in  eqs.  (1 ) ,  the  Q3D  formulation  imposes 
the  condition  =  0.  But  in  a  real  laminate,  the  condition  le^y  “  ° 

Ay 

cannot  be  imposed  except  at  the  ends  of  the  laminate,  where  eqs.  (1)  are 
not  exact.  Hence,  eqs.  (1)  make  the  approximation  that  -  0  through 

out  the  laminate. 

Equations  (1)  also  impose  the  condition  of  zero  curvature  about  the 

AviQ  T  A  K  =  -  ^  =  0.  This  condition  is  closely  approximated 

X  3x^ 

when  a  specimen  is  clamped  in  a  test  machine.  Hence,  imposing  =  0  m 

an  analysis  is  reasonable. 

Each  ply  was  idealized  as  a  homogeneous,  elastic  orthotropic  material 
with  the  following  properties: 


11  ^ 

134 

GPa 

(19.5  X 

10® 

psi) 

‘22  " 

E33 

=  10.2 

GPa 

X 

00 

. 

10® 

psi) 

^2  ^ 

^13 

=  G23 

=  5.r>2  GPa 

(0.8  X 

10®  1 

psi) 

1  ^ 

V 

= 

=  0.3 

12 

13 

23 

These  properties  are  representative  of  those  often  assumed  for  graphite/ 
epoxy.  The  subscripts  1,  2,  and  3  correspond  to  the  longitudinal, 

transverse,  and  thickness  directions,  respectively,  of  a  zero-degree  ply. 
Also,  the  uniform  axial  strain,  e^,  was  arbitrarily  set  equal  to  0.001 
throughout  the  study. 

ANALYSIS 

Two  types  of  analyses  were  used:  Q3D  finite  element  analysis  and  a 
laminate  theory  technique  for  calculating  total  strain-energy-release 
rate,  G^.  The  Q3D  analysis  was  used  to  obtain  free-edge  stress  distribu¬ 
tions  and  components  of  the  straln-energy-release  rates,  Gj,  Gjj,  and 
®III‘  laminate  theory  technique  was  used  because  it  is  inexpensive 

and  yields  an  interpretation  of  the  effect  of  laminate  thickness  on  G^. 


The  theoretical  basis  of  Q3D  finite  element  analysis  may  be  found  in 
reference  6,  so  only  the  application  of  the  Q3D  analysis  to  infinitely 
thick  laminates  and  the  finite  element  idealization  will  be  discussed. 

Infinitely  Thick  Laminates 

The  key  aspects  of  ar.al>_is  of  infinitely  thick  laminates  are; 

(a)  selection  of  the  repeating  unit  and  (b)  application  of  boundary 
conditions.  Consider  the  infinite  laminate  [ ( 45/0/-45/90)  ]  shown 

S  OD 

schematically  in  Fig.  2.  A  representative  repeating  unit,  (45/0/-45/90) 

s 

is  drawn  with  solid  lines.  Only  a  single  repeating  unit  need  be  analyzed 
If  appropriate  boundary  conditions  can  be  identified  for  the  upper  and 
lower  surfaces  of  the  unit.  Note  that  in  the  infinite  laminate  the  inter¬ 
face  between  adjacent  units  is  a  plane  of  symmetry.  Consequently,  if  we 
impose  the  condition  g  =  0  along  the  upper  and  lower  surfaces  of  a 
single  repeating  unit,  that  single  unit  behaves  as  though  it  is  within  an 
infinitely  thick  laminate.  In  the  finite  element  analysis  the  condition 
3y  =  0  along  a  line  is  implemented  by  constraining  all  nodes  along  that 
line  to  have  the  same  w  displacement  (but  with  the  net  z  force  =  0). 

In  all  cases  the  repeating  unit  is  a  symmetric  laminate;  hence,  in 
practice  only  half  of  the  repeating  unit  was  actually  modeled  if  there 
were  no  delaminations  or  the  delaminations  were  symmetrically  located 
about  the  midplane.  Otherwise,  the  entire  repeating  unit  was  analyzed. 

The  same  procedure  is  used  whether  or  not  the  laminate  has  a  delami¬ 
nation.  However,  note  that  placing  a  delamination  in  the  repeating  _nit 
causes  every  similar  interface  in  the  infinite  laminate  to  have  the  same 
delamination. 


Finite  Element  Idealizations 

Because  of  symmetries  in  the  problems  shown  in  Figure  1 ,  only  one- 
quarter  of  a  representative  x  =  constant  plane  was  modeled.  Figure  3(a) 
shows  a  typical  model  for  the  edge  stress  analysis  of  a  symmetric  8-ply 
laminate.  The  model  has  792  nodes  and  736  four-node  isoparametric  quadri¬ 
lateral  elements.  Models  for  other  laminate  thicknesses  had  the  same 
degree  of  mesh  refinement.  All  the  models  had  a  laminate  widti  to 

thickness  ratio,  2b/t,  of  at  least  15. 

For  both  finite  and  infinitely  thick  laminates  the  displacement  func¬ 
tions  U  and  V  were  prescribed  to  be  zero  along  the  line  y  =  0  and 
VJ  was  prescribed  to  be  zero  along  the  line  z  =  0.  For  the  infinitely 
thick  laminates  all  nodes  lying  on  the  line  z  =  4h  were  constrained  to 
have  the  same  displacement  W,  with  the  net  force  in  the  z-direction  on 
these  nodes  set  equal  to  zero. 

Figure  3(b)  shows  a  typical  model  for  a  16-ply  laminate  with  a 
delamination.  The  laminate  width  to  thickness  ratio,  2h/t,  was  at  least  5 
for  all  models.  The  delamination  length  to  ply  thickness  ratio,  a/h,  was 
at  least  20.  Around  the  delamination  tip,  square  elements  with  sides  of 
length  h/4  were  used,  strain-energy-release  rates  Gj,  Gjj,  and  Gjjj. 
were  calculated  using  the  virtual  crack  closure  method  (ref.  7). 

Classical  Laminate  Theory  Analysis 

The  laminate  theory  technique  for  calculating  total  strain-energy- 
release  rate,  G^,  is  bajed  on  dividing  a  laminate  into  two  types  of 
regions:  in  one  type  the  strain  energy  varies  linearly  with  delamination 

length,  and  in  the  other  type  the  strain  energy  is  independent  of  delami¬ 
nation  length  (ref.  8).  A  long  delamination  is  assumed  (i.e.,  a  delami- 
nation  length  of  more  than  twice  the  laminate  thickness). 
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Consider  the  laminate  in  Fig.  4,  which  is  divided  into  four  regions. 
Regions  I  and  III  are  away  from  geometric  discontinuities,  i.e.,  either 
the  crack  tip  or  the  free  edge.  Hence,  the  stress  distributions  in  these 
regions  are  defined  by  laminate  theory  and  are  independent  of  delamination 
length.  The  strain  energy  in  regions  I  and  III  vary  linearly  with  delami¬ 
nation  length  (since  the  volume  in  each  region  changes  linearly).  Regions 
II  and  IV  have  complex  stress  distributions  becau^  .  the  crack  tip  and 
free  edge,  respectively.  But  the  stress  distributions  are  independent  of 
delamination  length.  Also  the  volume  in  each  region  stays  constant. 

Hence,  the  strain  energy  in  regions  II  and  IV  is  independent  of  delami¬ 
nation  length. 

In  this  study  was  held  constant.  Hence,  is  given  by  the 

change  in  strain  energy  in  regions  I-IV.  But  as  previously  discussed, 
only  the  strain  energy  in  regions  I  and  III  vary  with  delamination  length. 
Reference  2  showed  that  this  change  in  strain  energy  could  be  expressed  as 

e/t 

“  2n  “  ®III^  <5) 


where  Ej  and  the  axial  Young's  moduli  of  regions  I  and  III 

respectively  and  n  is  the  number  of  delaminations. 

In  reference  2  the  axial  Young's  moduli  are  calculated  using  CLT 
implicitly  assuming  ^  ^  ^xy  example  the  axial 

modulus  Ej  would  be 


where  x^ ^  is  the  (1,1)  element  of  the  inverse  extensicnal  stiffness 
matrix  for  region  I,  t  is  the  thickness  of  the  laminate. 
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For  some  laminates  equation  (6)  is  sufficiently  accurate,  as  demonstrated 

in  reference  2.  However,  if  the  laminate  or  sublaminate  exhibits  strong 

extension-bending  coupling,  the  above  technique  requires  modification. 

This  modified  technique  still  uses  CLT  to  calculate  the  axial  modulus,  but 

the  calculation  is  performed  assuming  N  =  =  <  =  m  =  =  0  and 

y  xy  X  y  xy 

These  assumptions  are  consistent  with  those  made  in  the  Q3D 
finite  element  analysis.  The  effective  axial  modulus  E  is  then 


E  = 


(7) 


The  laminate  theory  analysis  was  also  used  to  analyze  infinitely 
thick  laminates  synthesized  by  stacking  "repeating  units."  The  procedure 
is  conceptually  the  sam^  as  used  in  the  Q3D  analysis;  the  difference  is  in 
the  implementation.  In  the  laminate  theory  analysis,  lines  of  symmetry 
between  repeating  urdts  are  accounted  for  by  imposing  zero  curvatures 

and  For  infinitely  thick  laminates,  the  thickness  t  is 

that  for  the  repeating  unit. 


RESULTS  AND  DISCUSSIONS 

First,  the  interlaminar  stresses  near  the  free  edge  of  undelaminated 
laminates  will  be  discussed.  Presumably  delamination  initiation  is 
governed  by  these  stresses.  The  magnitudes  of  these  stresses  near  the 
free  edge  are  well  known  to  depend  on  the  finite  element  mesh  refinement. 
In  fact,  the  stresses  are  singular  (mathematically)  at  the  intersections 
of  the  free  edge  and  the  ply  interfaces.  But  the  finite  values  calculated 
by  the  finite  element  analysis  do  estimate  the  intensity  of  the  stress 
field.  Hence,  the  calculated  stresses  do  qualitatively  show  the  effect  of 
laminate  thickness  on  interlaminar  stresses.  Secondly,  laminates  with 
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delaminations  are  considered.  Strain-energy-release  rates  vdll  be  used  to 
characterize  the  effect  of  laminate  thickness  on  interlaminar  stresses 
near  the  tip  of  a  delamination.  Finally,  a  few  experimental  results  from 
the  literature  will  be  interpreted  based  on  the  analytical  study. 

Pr^e-Edge  Interlaminar  Stresses 

Free-edge  interlaminar  stresses  depend  on  differences  in  the 
mechanical  properties  of  the  various  plies  in  a  laminate  and  the  stacking 
sequence.  For  thin  laminates  simple  equilibrium  arguments  can  be  used  to 
show  the  effect  of  stacking  sequence  on  o^.  For  example,  using  equilib¬ 
rium  arguments  the  midplane  for  a  (45/0/-45/90)g  laminate  is  pre¬ 

dicted  to  have  a  large  tensile  value  for  tensile  applied  load.  Similar 
arguments  predict  a  large  comprcsssive  at  the  midplane  for  the 

reversed  stacking  sequence,  (90,/-45/0/45)g.  But  the  infinitely  thick 
counterparts  for  both  of  these  laminates  is  [ (45/0/-45/90)g]„;  th.'t  is, 
l{45/0/-45/90)gJ„  is  the  same  as  I (90/-45/0/45 )gl«.  Obviously,  the  effect 
of  ” stacking  sequence”  becomes  clouded  for  thick  laminates. 

To  investigate  the  effect  of  laminate  thickness  on  free-edge  inter¬ 
laminar  stresses,  several  laminates  synthesized  by  stacking  ply  groups  of 
{45/0/-45/90)  were  analyzed.  The  laminates  considered  were: 

( 45 /O /-45 /90 ) g ,  ( 90/-45/0/45 ) g ,  I ( 45/0/-45/90 ) gloo#  ^  45/0/-45/90/90/-45/0/45 ) g , 
( 45 /O /-45/90/45/0/-45/90 ) g ,  and  [ ( 45/0/-45/90/45/0/-45/90 ) glo® •  First,  the 
effect  of  laminate  thickness  on  free-edge  stresses  for  interior  ply  groups 
was  examined.  Eight-ply  and  infinitely  thick  laminates  were  used  in  this 
first  phase.  Next,  two  16-ply  laminates  were  examined  to  determine  how 
laminate  thickness  affects  the  response  of  exterior  ply  groups  and  whether 
the  interior  ply  groups  of  a  1 6— ply  laminate  behave  as  though  they  were 
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inside  an  infinitely  thick  laminate.  Finally  the  effect  of  laminate 

thickness  on  free-edge  boundary  layer  width  was  examined. 

The  effect  of  laminate  thickness  on  the  free-edge  stresses  for 

interior  ply  groups  was  determined  by  analyzing  three  related  laminates; 

(45/0/-45/90)g,  (90/-45/0/45)g,  and  I (45/0/-45/90)3]„.  Note  that  because 

of  symmetries,  t  ( 45/0/-45/90)g]^„  .is  identical  to  [  (90/-45/0/45)^ 

Hence,  the  infinitely  thick  laminate  is  related  to  both  8-ply  laminates. 

Figure  5  shows  the  distribution  of  a  and  a  through  the 

z  xz 

thickness  of  the  specimen  at  the  free  edge.  These  stresses  reach  their 

maximxim  calculated  values  at  different  interfaces;  is  large  at  z  =  0 

and  h  and  is  large  at  z  =  2h  and  3h.  The  infinitely  thick 

counterpart,  [(45/0/-45/90)3]„  has  very  similar  stress  distributions  (see 

Fig.  6).  The  distribu’^ion  for  the  infinite  laminate  is  shifted  in 

the  compression  direction,  but  the  shape  of  the  curve  is  similar,  except 

in  the  45  degree  ply  where  the  infinite  laminate  has  significant 

compressive  o^.  m  contrasst,  the  distributions  for  the  two 

laminates  are  nearly  identical.  Figure  7  shows  the  o  distributions 

z 

from  Figures  5  and  6  along  with  the  a  distribution  for  the 

z 

(90/-45/0/45)g  laminate.  Note  that  the  distribution  for  the 

infinitely  thick  laminate  is  approximately  an  average  of  the  a 

z 

distributions  for  the  two  related  8-ply  laminates.  The  infinitely  thick 

laminate  has  the  detrimental  characteristics  of  both  thin  laminates  - 

has  both  large  tension  and  compression  o^.  But  the  maximum  tensile  and 

compressive  values  are  less  than  those  in  the  two  8-ply  laminates. 

The  response  of  interior  and  exterior  ply  groups  was  studied  using  two 

16-ply  laminates;  (45/0/-45/90/90/-45/0/45)g  and  (45/0/-45/90/45/0/-45/90) 

Figure  8  shows  the  through-thickness  a  and  o  distributions  for  the 

z  xz 
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(45/0/-45/90/90/-45/0/45)^  laminate.  Stress  distributions  for  (45/0/-45/90)g 

and  [(45/0/'-45/90)j^loo  superimposed  on  the  16-ply  results.  The 

and  o  distributions  for  the  interior  plies,  i.e.,  z  =  0  to  4h,  are 
xz 

very  similar  to  that  for  [ ( 45/0/-45/90) in  which  all  the  plies  are 
interior  plies.  But  from  z  =  4h  to  8h  the  and  distributions 

are  very  similar  to  that  for  (45/p/-45/90)g,  in  which  there  is  no  interior 
group  of  (45/0/-45/90)  plies.  Apparently  and  are  governed  by 

the  neighboring  plies.  Furthermore,  the  effect  of  the  top  surface  (i.e., 

2  =s  8h)  does  not  propagate  far  into  the  laminate.  This  interpretation  was 
verified  by  examining  another  16-ply  laminate,  {45/0/— 45/90/45/0/— 45/90 )g. 

In  the  region  4h  4  z  <  8h,  the  and  distributions  are  nearly 

the  same  as  that  for  a  (45/0/-45/90)g  laminate  (see  Fig.  9).  Hence, 
(45/0/-45/90)5,  (45/0/-45/90/90/-45/0/45)g,  and  (45/0/-45/90/45/0/-45/90) ^ 
laminates  all  have  approximately  the  same  o  and  a  distributions  in 
the  outer  four  plies,  even  though  the  remaining  plies  are  much  different 
for  each  laminate.  Near  z  =  4h  in  the  (45/0/-45/90/45/0/-45/90)g 
laminate  (Fig.  9),  the  stress  distributions  are  not  matched  by  those  for 
( (45/0/-45/90)g]^.  This  is  because  the  stacking  sequence  near  z  =  4h  in 
the  16-ply  laminate  is  different  from  that  found  anywhere  in  the  infinite 
laminate.  But  near  z  =  0,  the  stacking  sequence  is  the  same  as  that  in 
the  infinite  laminate;  concomitantly,  the  distributions  agree  very 

well. 

The  effect  of  laminate  thickness  on  the  free-edge  boundary  layer 
width  was  also  examined.  This  effect  v;as  qualitatively  assessed  by 
comparing  the  variation  of  in  the  width  direction  near  the  free  edge. 

Figure  10  shows  the  through -width  distributions  for  several  laminates 

along  90-90  interfaces. 
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The  boundary  layer  width  is  less  for  an  infinitely  thick  laminate 
than  for  the  corresponding  finite  thickness  laminate.  That  is,  the  width 
is  less  for  [{45/0/-45/90)3]^  th«n  for  (45/0/-45/90)3  (Fig.  10(a)). 

Hence,  the  boundary  layer  width  does  not  necessarily  increase  with 
laminate  thickness.  Increasing  the  number  of  plies  in  the  repeating  unit 
can  increase  the  boundary  layer  width,  as  evidenced  by  the  results  for 
[ (45/0/-45/90)g)g,  which  has  8  plies  in  the  repeating  unit,  and 
(45/0/-45/90/45/0/-45/90)3,  has  16  plies  in  the  repeating  unit 

(Fig.  10(b)).  (Of  course,  for  this  16-ply  laminate  there  is  only  a  single 
"repeating"  unit.)  By  increasing  the  number  of  plies  in  the  repeating 

unit,  a  laminate  becomes  less  homogeneous,  which  causes  a  larger  boundary 
layer  width. 

Strain  Energy  Release  Rates: 

As  stated  earlier  two  types  of  analyses  were  used  to  calculate  strain 
energy  release  rates:  a  laminate  theory  analysis  and  a  Q3D  finite  element 
analysis.  First  the  laminate  theory  analysis  will  be  used  to  show  the 
effect  of  laminate  thickness  on  total  strain  energy  release  rates.  Then 
the  Q3D  analysis  will  be  used  to  show  the  effect  of  laminate  thickness  on 
the  relative  magnitudes  of  Gj,  Gjj,  and  G^jj.  The  strain-energy 
release  rates  are  divided  by  h  in  the  figures,  since  the  strain-energy 
release  rates  vary  linearly  with  h. 

Classical  Laminate  Theory  Results 

Delamination  can  occur  at  a  single  interface  or  at  several  interfaces. 
The  possible  combinations  of  delaminations  are  too  numerous  to  examine  them 
all  individually.  Instead  two  patterns  of  delaminations  were  considered:  a 
single  de lamination  at  z  =  nh  and  delaminations  at  z  =  ±nh.  The  two 
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de laminations  at  z  =  ±nh  were  assumed  to  extend  simultaneously.  Recall 
that  for  infinitely  thick  laminates,  all  repeating  units  have 
delaminations  at  similar  interfaces. 

Delaminations  at  z  =  nh  and  z  =  ±  nh  in  (45/0/-45/90)g  and 
[ {45/0/-45/90)g]^  were  studied  first.  For  the  infinitely  thick  laminate, 
each  repeating  unit  ^which  consists  of  8  plies)  has  the  delaminations 
indicated.  Figure  11  shows  the  results  for  these  laminates. 

Specimen  thickness  has  virtually  no  effect  on  the  total  strain- 
energy-  release  rate,  G,^,  when  the  delaminations  are  at  z  -  ±  h,  ±2h,  or 
t3h.  But  for  single  delaminations  at  z  =  0,  h,  2h,  or  3h,  there  are 
large  differences;  for  the  infinitely  thick  laminate  is  always  much 

smaller  than  for  the  8-ply  laminate.  The  most  striking  case  is  at  z  =  0 
where  G,p  is  identically  zero  for  the  infinite  laminate,  but  is  very 
large  for  the  8-ply  case. 

These  trends  can  be  easily  explained  by  examining  the  basis  of  the 
CLT  analysis.  Changes  in  axial  modulus  are  what  contribute  to  G^.  A 
delamination  at  z  =  0  in  the  infinitely  thick  laminate  simply  divides 
the  original  quasi-isotropic  laminate  into  many  individual  quasi-isotropic 
laminates.  But  the  axial  modulus  of  a  quasi-isotropic  laminate  is  inde¬ 
pendent  of  thickness  or  stacking  sequence.  Hence,  G^  is  zero.  For 
z  =  h,  2h,  and  3h  the  differences  are  not  as  dramatic,  but  still  a 
de lamination  in  each  repeating  group  does  not  affect  the  axial  modulus  of 
a  thick  laminate  nearly  as  much  as  it  does  for  a  thin  laminate.  Hence, 

G^  is  smaller  for  the  thick  laminate. 

Recall  that  Figures  5  and  6  showed  that  for  an  undelaminated  speci¬ 
men,  0^  and  distributions  along  the  free  edge  were  nearly  the  same 

for  the  8-ply  cind  the  infinite-ply  laminates.  But  Figure  11  shows  that 
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is  identically  zero  for  a  delamination  at  z  =  0  in  an  infinitely 
thick  laminate,  and  is  large  for  a  delamination  at  z  =  0  in  an 

8-ply  laminate.  Since  is  a  measure  of  the  intensity  of  the  stress 

field  around  a  delamination  tip,  there  appears  at  first  to  be  a  contradic¬ 
tion.  This  apparent  contradiction  is  a  result  of  the  way  varies  with 

delamination  length,  a,  at  the  midplane  as  shown  in  Figure  12.  (These 
results  were  obtained  with  the  Q3D  finite  element  analysis.)  Note  that 
Gt  increases  with  "a”  initially,  as  would  be  expected  for  a  delamination 
initiating  in  a  stressed  material.  But  with  further  increase  in  delamina¬ 
tion  length,  decreases  to  zero  and  remains  zero  for  further  growth. 

The  laminate  theory  technique  always  calculates  the  "steady-state"  value, 
which  in  this  case  is  zero. 

These  observations  indicate  that  there  may  be  a  correlation  between 
free-edge  interlaminar  stresses  before  delamination  initiation  and  strain- 
energy— release  rate  right  after  de lamination  initiation.  But  after  the 
delamination  grows  away  from  the  free  edge,  there  is  not  necessarily  any 
correlation  between  the  original  free— edge  stresses  and  the  strain— energy- 
release  rate. 

Next,  (45/0/-45/90/45/0/-45/90)  and  [ ( 45/0/-45/90/45/0/-45/90)  ] 
laminates  were  examined.  The  repeating  group  in  these  laminates  consists 
of  16  plies.  As  in  the  cases  studied  earlier,  G,p  for  the  infinite  lami¬ 
nate  is  less  than  or  equal  to  G,p  for  the  thinner  laminate  (see  Fig.  13). 

The  last  phenomenon  studied  using  the  laminate  theory  analysis  was 
the  effect  of  delamination  location  on  G,p.  Only  a  single  delamination  is 
assumed  to  exist.  The  laminates  examined  were  [ ( 45/0/-45/90)g]^,  n  =  2, 

4,  8.  Figure  14  shows  the  results  for  these  16-,  32-,  and  64-ply  laminates. 
A  delamination  far  from  the  top  surface  of  the  laminate  affects  the  axial 


stiffness  very  little.  Hence,  the  corresponding  is  small.  For  thin 

laminates,  e.g.,  the  16-ply  laminate,  all  of  the  plies  can  be  considered 
to  be  near  the  top  surface.  Concomitantly,  for  the  16-ply  laminate  there 
is  no  distinct  "interior  region"  in  which  is  small.  But  for  the  32- 

and  64-ply  laminates  there  is  a  distinct  region  in  which  G^  is  small. 
Note  that  for  these  three  laminates,  the  interface  which  corresponds  to 
maximum  Gip  is  the  fifth  interface  from  the  top  surface  (i.e.,  at 
z  =  (4n-5)h).  Furthermore,  the  maximum  value  appears  to  be  asymptotically 
decreasing  to  a  constant  value  as  the  laminate  thickness  increases. 
Apparently  the  behavior  (in  terms  of  G^)  of  the  plies  near  the  surface  is 
nearly  independent  of  the  laminate  thickness.  Recall  that  the  free-edge 
stresses  for  the  plies  near  the  surface  discussed  earlier  were  also  nearly 
in(dependent  of  the  laminate  thickness. 

These  results  predict  that,  based  on  total-strain-energy-release 
rates,  de laminations  would  be  expected  to  grow  most  easily  near  the  top 
surface  rather  than  the  interior  of  a  thick  laminate.  Also,  outer  plies 
of  a  thick  laminate  should  be  less  prone  to  delaminate  than  the  outer 
plies  of  a  thin  laminate.  Furthermore,  since  the  thick  laminate  has  a 
smaller  percentage  of  outer  plies  than  a  thin  laminate,  the  percentage  of 
delaminated  interfaces  should  also  decrease  as  the  laminate  thickness 
increases. 

Finite  Element  Results 

The  CLT  analysis  calculates  only  the  total  strain-energy-release 
rates.  To  determine  the  relative  magnitudes  of  Gj,  and  Q3D 

finite  element  analysis  was  used. 

Eight-,  16- ,  and  32-ply  laminates  and  infinitely  thick  laminates  with 
8-,  16-,  and  32-ply  repeating  groups  were  analyzed.  The  finite  thickness 
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In  contrast  to  the  results  above,  the  infinitely  thick  laminates 
^ Table  3)  have  essentially  no  mode  I  component.  As  mentioned  earlier, 
^III  negligible.  Hence,  the  strain-energy-release  rate  is 

essentially  all  mode  II. 

Comparison  of  Analysis  and  Experiments 

The  analyses  presented  earlier  were  used  to  discuss  a  few  relevant 
experimental  results  in  the  literature. 

The  static  unnotched  failure  strains  of  [45/0/-45/90]g  and 
[45/0/-45/90] graphite/epoxy  laminates  (T300/5208)  have  been  measured 
as  0.854%  (ref.  3)  and  1.03%  (ref.  4),  respectively.  The  specimens  tested 
in  reference  3  and  4  were  manufactured  about  the  same  time  by  the  same 
vendor.  The  [45/0/-45/90] ^  laraiiiate  showed  extensive  delamination  growth 
before  final  failure.  In  contrast,  the  [45/0/-45/90] 2g  laminate  did  not 
show  any  delaminations  prior  to  final  failure,  even  though  the  final 
failure  strain  was  much  larger.  Free-edge  stresses  and  atrain-energy- 
release  rates  were  examined  to  determine  whether  one  would  expect  the 
[45/0/-45/90]g  laminate  to  delaminate  and  the  [45/0/-45/90] jg  laminate  not 
to  delaminate.  Figures  5  and  9  show  both  laminates  have  large 
interlaminar  normal  and  shear  stresses  at  the  free  edge.  In  fact,  the 
peak  values  of  these  stresses  are  nearly  the  same  for  the  two  laminates. 
Hence,  the  free-edge  stress  distributions  do  not  explain  why  the  8-ply 
laminate  is  more  prone  to  delaminate  than  the  16-ply  laminate. 

Strain-energy-release  rates  are  not  strictly  applicable  to  delamina¬ 
tion  initiation,  since  the  G's  are  zero  for  zero  delamination  length. 

But  the  G's  increase  quickly  to  constant  values  as  a  delamination  grows 
(ref.  2).  Hence,  if  small  initial  flaws  are  assumed  to  exist,  the  G's 
calculated  herein  are  appropriate  for  predicting  delamination  growth. 
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The  strain-energy-release  rates  for  the  two  laminates  are  quite 
different.  Figures  11  and  13  show  that  for  the  same  strain  level,  G,j,  is 
much  larger  for  a  delamination  in  a  t45/0/-45/90] ^  laminate  than  for  a 
[45/0/-45/90] 2g  laminate,  for  both  single  and  two  symmetric  delaminations. 
Based  on  one  would  predict  a  single  delamination  to  grow  in  the  8-ply 

laminate  along  z  =  h. 

In  the  tests  (ref.  3)  a  single  delamination  did  form  in  the  8-ply 
laminate,  but  it  wandered  back  and  forth  from  one  -45/90  (i.e.  z  =  h) 
interface  across  the  90  deg  plies  to  the  other  —45/90  interface  (see 
Fig.  15).  This  pattern  of  growth  cannot  be  rigorously  analyzed  using 
either  the  Q3D  or  the  laminate  theory  analysis,  since  the  problem  is  fully 
three-dimensional.  The  actual  G,ji  was  assumed  to  fall  between  G,j, 
calculated  for  a  delamination  at  z  =  0  and  G,j,  for  a  delamination  at 
z  =  h,  since  the  actual  delamination  path  is  bounded  by  z  ==  ±h.  Since 
G^  for  a  delamination  at  z  -  0  is  larger  than  the  maximum  Gr^  for  the 
16-ply  laminate,  one  would  still  predict  the  8-ply  laminate  to  delaminate 

at  a  lower  strain  than  the  16-ply  laminate. 

The  possibility  that  delamination  growth  is  governed  by  the  magnitude 
of  the  mode  I  component  Gj  was  also  considered.  Except  for  midplane 
delamination,  for  which  =  G^,  the  Q3D  analysis  is  required  to 
determine  the  mode  I  component.  Q3D  analysis  was  performed  for  delamina¬ 
tions  at  z  =  h  and  z  =  ±2h  for  both  laminates.  The  16-ply  laminate 
was  also  analyzed  for  a  delamination  at  z  =  5h.  These  interfaces  were 
chosen  because  G^  was  large  at  these  interfaces.  Table  4  shows  the 
results.  Based  on  one  would  predict  the  8-ply  laminate  to  delaminate 

at  a  lower  strain  level  than  the  16-ply  laminate.  The  delamination  in  the 
8-ply  laminate  should  grow  along  the  midplane  (i.e.  z  =  0).  For  a 


delamination  in  the  8-ply  laminate  along  z  =  h  (where  is  largest), 

Gj  is  still  much  larger  than  that  for  any  of  the  16— ply  cases#  Assuming 
the  actual  Gj  is  bounded  by  the  results  for  delaminations  at  z  =  0 
and  z  =  h,  one  would  predict  (correctly)  that  the  8-ply  laminate  should 
delaminate  at  a  lower  strain  level  than  the  16— ply  laminate# 

The  critical  value  of  Gj  for  delamination  growth  (G^^)  in  T300/5208 
is  approximately  100  J/m^  (ref#  9)#  Assuming  that  delamination  growth 
occurs  when  Gj^  is  reached,  the  previous  analytical  results  can  be  used 
to  predict  when  delamination  should  occur  for  each  laminate# 


For  the  8-ply  laminate  the  average  value  of  Gj  for  delaminations 
at  z  =  0  and  z  =  h  was  used  (i.e#  G^/h  for  e  =  0#001  is 

A  O 

(20.62  +  13.48)/2  =  17.05  J/m^).  The  ply  thickness  h  was  assumed  to 

be  1,397  X  io“^  m#  The  strain-energy-release  rate  varies  as  e^# 

o 

Hence , 


^CR\ 

17.05  kJ/m^ 


100  J/m 


1  .397  X  10“\ 


(8) 


where  is  the  strain  at  which  delamination  should  occur#  Solving 

equation  (8),  one  obtains  =  0.65%.  In  reference  3  is  not 

o  o 

reported.  But  the  stress-strain  curve  is  reported  and  becomes  nonlinear 
(probably  due  to  delamination)  between  0.6%  to  0.7%  strain.  Hence,  the 
agreement  between  analysis  and  experiments  is  excellent. 

For  the  16-ply  laminate  the  largest  Gj  occurs  for  delamination  at 
z  =  0  (see  Table  4).  The  magnitude  of  G^/h  for  e  =  0.001  is 
9.79  kJ/m^.  Proceeding  as  for  the  8-ply  laminate,  is  predicted  to 

be  0.86%#  In  the  actual  test  the  specimen  failed  at  a  strain  of  1% 
without  delaminating.  Of  course,  whether  delamination  growth  was  imminent 
when  failure  occurred  is  not  known. 
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Compression-compression  fatigue  test  results  for  t90/45/0/-45] G/E 
laminates  have  been  reported  in  reference  10.  In  these  tests  the  first 
de lamination  was  observed  to  initiate  near  an  outside  ply.  Subsequent 
delaminations  formed  farther  in  the  interior  of  the  laminate  (i.e.,  closer 
to  the  midplane).  This  behavior  is  consistent  wich  the  analytical  results 
discussed  earlier;  i.e.,  the  total  strain-energy-release  rate  and  the 
percentage  of  Gj  tend  to  be  larger  for  de lamination  growth  near  the 
outside  plies. 


CONCLUSIONS 

Analyses  were  conducted  to  determine  the  effect  of  laminate  thickness 
on  the  interlaminar  stresses  in  rectangular  laminates  under  uniform  axial- 
strain.  Laminate  thicknesses  ranged  from  8-ply  to  infinitely  thick. 

All  the  laminates  were  quasi-isotropic.  They  were  synthesized  from 
(45/0/-45/90)  ply  groups.  For  laminates  without  delamination,  free-edge 
stresses  were  examined;  for  laminates  with  delamination,  strain-energy- 
release  rates  were  used  to  characterize  the  interlaminar  stresses  at  the 
delamination  tip.  Quasi-three-dimensional  finite  element  analysis  and 
classical  laminate  theory  analysis  were  used  in  the  study.  Based  on  the 
results  from  these  analyzes,  the  following  conclusions  were  made; 

1 •  If  thick  laminates  are  synthesized  by  stacking  repeating  units,  then 

(a)  the  interlaminar  normal  stress  distribution  in  the  outer  group  of 
plies  is  insensitive  to  the  total  laminate  thickness, 

(b)  the  interlaminar  normal  stress  distribution  for  the  interior 
group  of  plies  is  nearly  the  same  as  that  for  an  infinitely  thick 
laminate . 

(c)  The  interlaminar  shear  stress  distribution  is  very  insensitive  to 
the  total  laminate  thickness. 


2/ 


.  ••  /  \  1 


The  free-edge  boundary  layer  width  is  related  to  the  thickness  of  the 
repeating  unit,  not  the  total  laminate  thickness. 

Free-edge  interlaminar  stresses  can  be  used  only  for  locating 
delaraination  initiation.  After  a  delamination  grows  away  from  the 
free  edge,  the  original  free-edge  stresses  provide  no  information  on 
the  intensity  of  the  stress  field  at  the  delamination  tip. 
Delaminations  near  the  top  and  bottom  surfaces  of  a  thick  laminate 
have  larger  total  strain-energy-release  rates,  G^,  and  mode  I-to- 
total  ratios  than  delaminations  deep  in  the  interior  of  the 

laminate.  Therefore,  delaminations  can  be  expected  to  grow  most 
easily  near  the  surface  of  a  laminate.  This  expectation  is  consistent 
with  experimental  results  in  the  literature.  Also,  near  surface 
delaminations  in  thin  laminates  tend  to  have  larger  strain-energy 
release  rates  than  corresponding  near  surface  delaminations  in  thick 
laminates.  Therefore,  near  surface  delaminations  can  be  expected  to 
grow  more  easily  in  a  thin  laminate  than  in  a  thick  laminate. 
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LAMINATE* 


8-PLY 


1 45/0/-45/90/90/-45/0/45 ) , 


21  .77 


1 6-PLY 


[ 45/0/-45/90/90/-45/0/45/ ( 45/0/-45/90 )  ] 


11.24 


[ 45/0/-45/90/90/-45/0/45/ ( 45/0/-45/90 )  ] 


19.46 


32-PLY 


[45/0/-45/90/90/-45/0/45/{45/0/-45/90)  3] 


9.143 


T 

[ 45/0/-45/90/90/-45/0/45/ ( 45/0/-45/90)  3 ] 


16.29 


r 

[ {45/0/-45/90)g/45/0/-45/90/90/-45/0/45/(45/0/-45/90)  _] 


7.433 


[  (45/0/-45/90)g/45/0/-45/90/90/-45/0/45/(45/0/-45/90)  -] 


6.612 


INDICATES  DELAMINATION  LOCATION 


TABLE  1 .-  Strain-energy  release  rates  for  various  laminates 
vxS.  a  single  delamination. 


LAMINATE* 

Gj/Gt 

G^/h,  kJ/m^ 

8-PLY 

+ 

(45/0/-45/90]g 

.62 

12.82 

1 6-PLY 

+ 

[45/0/-45/90/90/-45/0/45)g 

.55 

9.474 

[45/0/-45/90/90/-45/0/45 ] g 

.22 

21  .93 

32-PLY 

4 

[ 45/0/-45/90/90/-45/0/45/ ( 45/0/-45/90 )  ] 

s 

.57 

8.468 

4 

[ 45/0/-45/90/90/-45/0/45/ ( 45/0/-45/90 )  ] 

®  s 

.452 

15.71 

[ (45/0/-45/90)g/45/0/-45/90/90/-45/0/45]g 

.01 

7.081 

4 

[ (45/0/-45/90)g/45/0/-45/90/90/-45/0/45]g 

.03 

11  .34 

*  "4'”  INDICATES  DELAMINATION  LOCATION 


TABLE  2,-  Strain-energy  release  rates  for  various  laminates 
with  two  symmetrically  located  delaminations. 
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LAMINATE 


8-Pt.Y  REPEATING  UNIT 


[(45/0/-45/90)g] 


16-PLY  REPEATING  UNIT 


T 

[ (45/0/-45/90/90/-45/0/45)3] 


T 

[ (45/0/-45/90/90/-45/0/45)g] 


32-PLY  REPEATING  UNIT 


[( (45/0/-45/90/90/-45/0/45/(45/0/-45/90)g)  ] 

®  0 

[ ( 45 /0/-45/90/90/-45/0/45/ ( 45/0/-45/90 ) g)  ] 


T 

[((45/0/-45/90)g/45/0/-45/90/90/-45/0/45)3]^ 


T 

[[ {45/0/-45/90)g/45/0/-45/90/90/-45/0/45)g]^ 


.057 


.0001 


.0001 


Gf^/hf  kJ/in 


.0005 


.03 


.03 


0.0 


12.95 


9.522 


9.522 


8.626 


3.337 


3.337 


8.626 


*  "4-”  INDICATES  DELAMINATION  LOCATION 


TABLE  3.-  Strain-energy  release  rates  for  various 
infinitely  thick  laminates. 
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Fig,  2 — Schematic  of  [ (45/0/-45/90)  ]  laminate. 
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b)  TYPICAL  EDGE  DELAMINATION  MODEL  443 
NODES,  384  ELEMENTS 


Fig.  3 — Typical  finite  element  models. 


Fig,  4  Laminate  divided  into  four  subregions,  I  -  IV 


ORSGiaAL  PAGi-  IS 
OF  POOR  QlIALifY 


Fig.  7 — Distribution  of  at  free  edge  for  [90/-45/0/45]^,  [ (45/0/-45/90)^] 

and  [45/0/-45/90]  laminates,  (e^  =  0.001) 


Fig.  8 — Comparison  of  and  distributions  at  free  edge  for  [45/0/-45/90/ 

90/--45/0/45]g  laminate  with  those  for  [45/0/-45/90]^  and  [ (45/0/-45/90)^ j 
laminates,  (e  «  0,001) 


/ 
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- 

•/ 

y!  r  - - 


ORiGif'Ai.  i  v;;  ,.  ,r: 
OF  POOR  QUkLlTi 


(45/  0/  -45/  90/  90/  -45/  0/  45)^  LAM  I  NATE 


Fig.  8 — Concluded. 
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]  ■  \ 


OF  POOR  Q 


m 


(45/  0/-45/90) 


(45/0/-45/90/45/  0/-45/  90). 


0  |[(45/  0/ -45/ 90) 


la 


a  MPa 

a)  0  DISTRIBUTIONS 


Fig.  9-Comparisou  of  0^  and  0^^  distributions  at  free  edge  for  a  [45/0/-45/90/ 
45/0/-45/90]g  laminate  with  those  for  [45/0/-45/90]g  and  [ (45/0/-45/90)g]^ 
laminates.  =  0.001) 


ORIGINAL  Pr  iP 
OF  PCCR  QL.'L'.r; 


Fig.  9 — Concluded. 
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a)  8-PLY  AND  INFINITELY  THICK  LAMINATES 


Fig.  10 — Through-width  distribution  of  along  SO-90  interface. 


DC  B  A 

II  II 

(45/0/ -45/90/90/ -45/0/45) 


DELAMI  NATION 


LOCATION 


ORiGf.-iAl  Pf 
OF  POOR  QU 


0/  -45/90) 


(45/  01  -45/90) 


ZERO 


G^/h,  Kj/m 

a)  SINGLE  DELAMI  NAT  I  ON 


Fig.  ll~Total  strain-energy-release  rate  for  various  delanination  locations  in 
[45/0/-45/90]  and  t(45/0/-45/90)g]^  laminates,  (c^  =  0.001) 
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OF  POOR  QUALITY 


C  B  A  ABC 

I  i  i  \  \  \ 

(45/0/ -45/90/90/ -45/0/45) 


b)  TWO  SYMMETRICALLY  LOCATED  DELAM I  NATIONS 


Fig,  11 — Concluded. 
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GF  EDCB  A  ABCDE  FG 

i  +  ll  f  i  WIf 

(45/0/ -45/90/45/0/ -45/90/90/  45/0/45/90/ -45/0/45) 


Cn:Ci:l:.L  v. 

OF  POOR  QUALITY 


DELAM I NATI ON  ^ 0/  -45/ 90)J ^ 


0  5  10  15  20  25 

G^/h,.kJ/m^ 

b)  TWO  SYIVIIVIETRICALLY  LOCATED  DELAM  I  NATIONS 


Fig.  13 — Concluded. 


45 


Of--  POOR  QUALITY 


Fig.  14 — Total  strain-energy-release  rate  vs.  distance  of  a  single  delamination 
from  top  surface,  (t/2  -  z)/h.  (e^  •»  0.001) 
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